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INTRODUCTION 


Recently, Ip and Goertz (1983), hereafter referred to as 
paper I, have suggested that the observed local time asymmetry of 
the UV emission from the Io torus (Sandel and Shemansky, 1982) could 
be accounted for if there exists a dawn to dusk electric field at 
the orbit of Io. The magnitude of this field was estimated as 
1.5 mV/m. Even though this field may seem small it is not easy to 
understand that such a field can exist at all. As mentioned in 
paper I the usual solar wind induced polar cap electric field (if 
it exists) would correspond to a dusk to dawn electric field inside 
the Jovian magnetosphere, i.e., opposite to the solar wind induced 
convection field at the earth. This is simply due to the fact that 
the direction of the magnetic dipole at Jupiter is opposite to that 
at the earth. On the other hand, it seems very difficult to imagine 
a local time asymmetry without a solar influence (solar wind induced 
or otherwise). A possible mechanism for maintaining a dawn to dusk 
electric field is described in this paper. In essence this electric 
field is due to the fact that in the outer magnetosphere the 
centrifugal current caused by the rotation of magneto spheric plasma 
flows across a density gradient and is, hence, not divergence free. 
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A Birkeland current system results with currents into the ionosphere 
at dawn and away from the ionosphere at dusk. The current closure 
in the ionosphere sets up a dawn to dusk resistive electric field 
which is mapped out into the magnetosphere. How big can such a 
field be? 

It has been suggested by Krimigis et al. (1981) that the 
presumably Io injected plasma is confined to a disk whose equatorial 
cross section is not circular but resembles an ellipse. This can be 
understood when one considers the solar wind induced distortion of 
the Jovian magnetic field. On the front side, field lines are 
compressed and they are extended on the night side. The corotation 
induced electric field is then also asymmetric (Goertz, 1978) . The 
equipotential contours come closer to the planet on the day side. 
Since the flux tube content is constant along the drift path (Hill 
et al., 1981), constant flux tube contours must also be quasi- 
elliptical. 

Krimigis et al. (1981) and others before them (Michel and 
Sturrock, 197^5 Kennel and Coroniti, 1975> 1977 ; Hill et al., 197^-5 
Hill and Michel, 1976) have suggested that as the plasma convects 
into the night side field lines extending beyond a certain critical 
distance may break open (Vasyliunas, 1983) and a rotationally 
driven antisunward directed plasma wind should flow in the tail. 

In paper I we have suggested a direct connection between this wind 
and the dawn to dusk electric field. However, it is not clear that 
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the wind occurs on "closed" field lines or "reconnected" open field 
lines as suggested by Vasyliunas ( 1983 ), in which case the electric 
field inside the boundary between the rotational flow and the wind 
would be dusk to dawn corresponding to a sunward flow. However , it 
is plausible that flux tubes which cannot enforce rotation in the 
tail and from which plasma is lost will have a significantly reduced 
flux tube content. This is implied in plate 1 of Krimigis et al. 
( 1982 ). This flux tube content gradient will also exist in the front 
side if the filling by outward radial diffusion is not fast enough. 

This model suggests that there exists a fairly well defined 
boundary between the full and empty flux tubes. The equatorial trace 
of the "last" full flux tube is elliptical, as shown in Figure 1. 

In order to simplify the problem we will treat this boundary as 
infinitely thin. Inside the boundary the full flux tubes move on 
equipotential contours which must be calculated self -consistently. 
Outside, the empty flux tubes move on equipotentials which are given 
by the solar wind induced magnetic field distortion and the corre- 
sponding electric field. We also make the simplifying assumption 

that near this boundary the ionosphere cannot exert a torque on the 

2 

plasma, i.e., that ft* l/r . Hill (1979) has shown that this is a 

good approximation beyond a distance L where L » 20 is a critical 

o J o 

L-shell parameter related to the mass injection from the Io torus. 

There is, however, an azimuthal current j at the boundary due to the 
difference of the polarization drifts of low energy ions and electrons. 
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This current has a component towards the dense regions on the dusk 
side and towards the empty regions on the dawn side (see Figure l) . 
There is also a current due to the pressure gradient. The component 
of the pressure gradient, normal to the boundary, drives a current 
along the boundary which does not require a field -aligned current. 
There is, however, a component of the pressure gradient along the 
boundary because the plasma adiabatically cools as it convects into 
the decreasing magnetic field of the tail. This component drives a 
current normal to the boundary which has the same sign as the normal 
component of the centrifugal current. However, the plasma flow 
velocity (rotation) is most likely supersonic and the pressure 
gradient current is small compared to the centrifugal current. We 
shall neglect it and treat only the cold plasma case. 

The polarization current has to be fed by a Birkeland current 
system with a field-aligned current out of the ionosphere and into 
the equatorial plane on the dusk side and in the opposite direction 
on the dawn side. The closure of this current in the ionosphere 
requires a dawn to dusk Pederson current and hence a dawn to dusk 
electric field. This self-consistent electric field distorts the 
drift orbits of the full flux tubes in a way to be calculated below 
(Fejer, 1964; Vasyliunas, 1970). The expected result is that in 
the night side, flux tubes move radially outwards and on the day 
side they move inwards. When this motion is superimposed on the 
rotational motion a dawnward displacement of the drift path results. 
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We would then expect the .plasma to be colder at dawn than at dusk, 
as described in paper I. 

Before we proceed with the calculations, we would like to 
discuss some observational evidence which demonstrates the existence 
of the boundary postulated above. Figure 2 shows data obtained by 
the LECP instrument during the inbound path of Voyager 2 (Krimigis 
et al., 1981). Between the last magnetopause crossing at a distance 
of 63 Rj and about 45 Rj, the LECP density, i.e., the density of hot 
ions (T ~ 30 - 40 keV), balances practically all the electron 
density obtained from the plasma wave (PWS) data (Gurnett et al., 
1981 ) . The density of the lower energy ions , determined by the PLS 
instrument (Bridge et al., 1979a, b; Belcher et al. , 1980; and 
McNutt et al., 1981), is quite low in this region. Inside of about 
45 R t the plasma density in the current sheet is an order of magni- 
tude higher. Such a region of low density but high temperature 
plasma between the magnetopause and the fully developed (Gledhill) 
plasma disc can also be observed on the Voyager 1 and 2 outbound 
paths (Figure 3), as well as in the Saturnian magnetosphere (Goertz, 
1983) • It may be a characteristic feature of rapidly rotating 
magnetospheres in general (Goertz, 1983) . The observation on the 
night side of a planetary wind by Krimigis et al. (1981) is obviously 
strong evidence for the existence of the boundary. In the region 
it was observed by Voyager 2 (between l40 and l60 Rj) the PSW 
total electron density was reduced by at least an order of magnitude. 
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as compared to the density in the nightside Gledhill disc (see 
Figure 3). 
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I. GOVERNING EQUATIONS 

In the magnetosphere the plasma is confined to a thin 

equatorial Gledhill disc across which neither the magnetic field 

(B = - B (r)e ) nor the plasma velocity changes significantly. The 
o z 

mass density per unit magnetic flux T], 

W?** - <D 

is constant along the drift path (Hill et al., 1981 ) 

V E • V71 = ( B -* - 2 — ) • VT1 = 0 , (2) 

-4 

where the electric field E is expressed in terms of the electrostatic 
potential § . To first order the drift velocity V^ is in the 
azimuthal direction, i.e., 

E - - V§ = - flRj sin 6 X B - VU , (3) 

where U is a small "resistive" potential to be calculated below. 
Associated with the plasma drift is a polarization current which, 
when integrated across the Gledhill disc, is given by 
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J X = H b X (V E * V V E ) 


( 4 ) 


Note that we have assumed the plasma to he cold and consequently 

-4 

neglected the pressure gradient current. Since V • J = 0 the 

■4 ( 

divergence of Jj. must be cancelled by a field-aligned current j 
(Vasyliunas , 197 0 ) , 



This current maps into a field -aligned current into the ionosphere 


. i 


j , which is given by 




( 6 ) 


where B 1 is the magnetic field strength in the ionosphere. Current 
continuity in the ionosphere requires (superscript i refers to the 
ionosphere) 


sin k 1 = V X 1 [S 1 (E + v n X B 1 )] 


( 7 ) 


The magnetic field inclination is h 1 and F 1 is the height integrated 
conductivity tensor. The neutral gas velocity is v^ = fijBj sin B 
In the equator the drift velocity is 


<(U 


w | « Hr 


(8) 


-» 



fir e + w 
V 


> 


with 


w = VU X e /B , 
z' o 


and hence neglecting terras which are quadratic in w. 


J = T] fi 2 r e + T]{w (2fi + r|% 

cp L r dr 


5W 


fi"r^ } e 
ocp r 


3w 

+ T]{2fiw cp - 


(9) 


The first term is the well-known centrifugal current, the second, 
term represents a radial current which exerts a torque on the 
radially drifting plasma, and the third term is a small correction 
to the first, which we shall neglect. 

The radial current has been treated ty Hill (1979) who has 
shown how it is related to the mass flux injected into the magneto- 
sphere from the Io torus. We will use his solution of the equation 


• • 

V x [Zr(fi T - fi)R, sin 9 

= J J 


X B 1 ] = 


sin h 1 Ib 1 1 1 5 / 


( 10 ) 


which predicted quite accurately the observed (Bridge et al., 1979) 
decrease of fi with radial distance. Hill has shown that beyond a 


certain distance, r Q = L R , the ionosphere cannot exert enough 

torque to maintain corotation. In fact, beyond r Q a very good 

2 

approximation is given by fl « l/r in which case » 0. Thus we 
write for the outer magnetosphere 


V i 1 (E i V 1 U) 


sin 1 B^~ | 2 

2 B “ 
o 



( 11 ) 


This equation is similar to equation l4 of Hill et al. (1981) who 
have used the same derivation albeit in a corotating frame of 
reference. In their work, cp is the longitude (e.g.. System III) in 
a frame rotating with Jupiter. In our case, cp is the longitude in 
an inertial, sun -fixed coordinate system, i.e., local time. 

Since V^, • 7T) = 0, the local time asymmetry of the flux tube 
content is related to the radial gradient 




- w |"(T]) 
r dr 


( 12 ) 


To solve equation 11 we need to map its right hand side to the 
ionosphere. To do that, we note that the potential is constant along 
a field line which we label by its Euler potential ct. In the equator 
the magnetic field is purely in the z -direction and given by 


r 


r or cp 



(13) 


11 


In the ionosphere we approximate the magnetic field by a pure dipole 
field with the dipole axis along the z-axis (rotation axis). We 
then have 


B r dr = - 2R T ^ sin 0 cos 9 B d9 
o J s 


(14) 


where B is the surface dipole magnetic field strength at the equa- 
s 

tor. Using the well-known equations for B 1 and sin k 1 for a dipole 
field and assuming, for simplicity, that £ is everywhere the same, 
we finally arrive at the equation 


1 9 / - n aux 1 

“7 a tT (sin 9 rr-j + 
sin 0 39 38 


. sfu 


o ¥ + g(B) 
j „2 a ^2 3cp 


sin^ 9 3 cp £ 


= 0 


(15) 


g( 9 ) = 


1 ~ 0 is (T| ) 


2E B sin' 0 
P s 


This equation must be solved subject to the boundary conditions to 
be discussed in the next section. 
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II. BOUNDARY CONDITIONS AND SOLUTIONS 


Like other authors before (Hill, 1979; Hill et al., 1981) , 
we will assume that within the plasmadisc plasma transport occurs 
via an adiabatic interchange diffusion which maintains a constant 
flux tube content T|. Outside the boundary, T) is much smaller and, 
for simplicity, we will set it equal to zero. The boundary field 

-K- ^ 

lines intersect the ionosphere at a latitude 0 . The value of 9 
is not well known. We will treat this as a free parameter. Inte- 
grating equation 15 with respect to 9 yields the first boundary 
condition 


[au/39 + 


cos 9 
sin 9 


Q 71 dU/59 j 

2Z B sin 7 9 
P s 


9=0 


* 


= 0 


(16) 


where [ ] indicates the difference between the value of the bracketed 
quantity inside and outside the boundary. In addition, we require 
the potential to be continuous across the boundary 

[u] e=0* = 0 • (l7) 

For 0 < 9 we have T] « 0 and it seems plausible that the ionosphere 
can enforce rigid corotation on the empty flux tubes at lower 
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colatitudes as 8 -» 0. In fact, the apparent disagreement between 
the LECP observations of nearly exact corotation in the low density 
region (Krimigis et al., 1981 and Carbary et al. , 1981) , and the PLS 
observations of significantly reduced rotation rates in the high 
density regions (Bridge et al., 1979a, b; McNutt et al., 1981) , indi- 
cates precisely this. Thus, in the low density region U -+ 0 as 
9 0. 

Since g(0) = 0 here the potential U must be of the form 

U° = ban 75 (2A° cos cp + 2B° sin cp) 

+ S ban m (0/2) (2 m A° cos rrsp + 2 m B° sin mcp) , (l8) 

m>l 

we will regard only the first harmonic terms (m = l) because we do 
not expect the solar wind to cause higher harmonic distortions of the 
magnetic and electric field. We further note that the electrostatic 

potential $ can be expressed as a function of the Euler potential O', 
say I = J {2(ct , )dQ! / . In that case, a is a constant along the drift 

contour and the form of Q : given by Stern (1967) together with equa- 
tion (3) can be used to derive approximate values for the coeffi- 
cients A° and B°. Since the noon -midnight meridian is the symmetry 
axis for the solar wind induced perturbations, we obtain B° = 0. The 
remaining v alu e of A° depends on the magnitude of the solar wind 
induced compression of the magnetosphere, i.e., the solar wind 
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stand-off distance, r = L R T , or equivalently the colatitude 

S S J 

-1/2 

sin 0 » L 'of the polar cap. Using the relation of Stern 


(1967) based on the magnetic field model of Mead (1964), we obtain 


A° “■ 0.27 R T 2 B fi T sin 0 6 
1 J o J s 


This corresponds to a very small poleward displacement of a field 

line by A8 ~ 16 L ^ degrees as it rotates from noon to midnight, 
s 

Inside the boundary the solution for U is of the form 
U 1 = tan(0/2) ^Ag 1 cos cp + 2 B 2 1 sin cp) 

+ tan" 1 (0/2)(2C 2 1 cos cp + 2D2 1 sin cp) . (20) 


Using the boundary conditions (16) and (17) we can solve for the 
potential U . A particularly simple form results when we approximate 
tan 0/2 by 0/2, which should be good for high latitudes. In that 



A sketch of this potential using an unrelatis tic ally low value for 7) 
is shown in Figure 4. Near Io's colatitude (0 ^ q = 26°) and with a 
realistic value for T|, the contours given by equation (21) are 


15 


indistinguishable on the scale of the drawing from those which 
result when the full function (tan 9/2), rather than its approxima- 
tion, is used. The compression of the equipotential contours between 
0 and 8^ o is very large for realistic values of 7] and it is diffi- 
cult to show individual contours clearly. Two points are important, 
(i) The electric field inside the boundary between full and empty 
flux tubes may be significantly larger than the original solar wind 
induced field. This depends on the value of the quantity 


x = 


Tin 


), . 7 2 B 

4 sm 0 p s 


*2 *2 


Whereas the potential difference (between noon and midnight) at 0 
is A° 0*, it is A° 0 ijl "+ x 2 at the colatitude 0. As we will see 
this may become quite large because x q = 65 . In Figure 4 we have 
used x q = 4 to make the figure clearer. 

(ii) The direction of the electric field changes as one goes to 
lower latitudes approaching a dawn to dusk direction, as 0 becomes 
larger. The angle of rotation increases with increasing x. 

At the inner edge of the torus we have again a steep radial 
gradient of T) (Goertz and Ip, 1982) and the electric field will 
change direction and magnitude. Whereas in the dense regions of the 
torus the plasma flow has an antisunward component, it is mainly 
sunward inside the inner edge of the torus. We have shown 


15 '. 


schematically the equipotential contours inside the inner torus 
boundary by dashed lines in Figure 4 . (These lines are not calcu- 
lated exactly.) 
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III . DISCUSSION 


Whereas the general result of the previous section is quite 

encouraging, it is not clear that a 1.5 mV/m field results at the 

orbit of Io, because 9 is probably quite small. However, Connerny 

s 

et al. (1981) estimated the magnetic latitude 0 of the auroral field 

s 

lines to be about ® s ^ 18° which would correspond to a dipole L 

shell parameter of only 10. Previous estimates, not including the 

current sheet, put 0 ^ 10° corresponding to L 33. We will use 

s s 

the lower value for 0 to he on the conservative side. A value for 

s 

-)f 

0 is equally difficult to estimate. Clearly, © s < 0 < 0 . One 

may argue that the value for 9 , estimated by Connerny et al. (1981), 

s 

really corresponds to 9 because, in essence, Connerny et al. have 
determined the footprint in the ionosphere of the field lines bound- 
ing the current sheet which is filled with dense plasma. We thus 
take 9 = 18 °. At this distance, the plasma rotation rate is reduced 

to something like 0.5 of its rigid value (McNutt et al., 1981) . The 
flux tube content in the outer magnetosphere is at least 
7) = 2 X 10 ^ kg/W (Hill et al. , 1981) but it may be as large as 
T) = 4 X 10 ^ kg/W, the value observed in the Io torus (Hill et al., 
1981) . In fact, using new results from Bagenal (private communica- 
tion) , one finds Tj = 8 X 10 ^ kg/W in the torus . Thus , ’Tj = 10 ^ kg/W 
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seems a reasonable value. The value for E is also not well known 

P 

(see Hill et al., 1981 for a discussion), but E^ = O.25 mho seems a 
reasonable number, even though values 10 times larger or smaller may 
also be realistic. 

With these numbers we get A° = 2500 V, x =“ 65 and at the 

1 o 

latitude of Io's flux tube, x =“ 35* Thus, the maximum of the poten- 
tial U occurs at a local time of 17:58, i.e., is almost exactly at 

dawn. The magnitude of the electric field at Io's orbit, however, 

_2 

is only 2 x 10 mV/m, i.e, a factor of 7 smaller than the field 
estimated in paper I. However, we do not regard this discrepancy as 
very serious because of the following two reasons; 

1. The value of x may be too small. There are at least two 

reasons for this. The value for T| used (10 ^ kg/W) may quite 
possibly be too low by a factor of eight. In the Io torus, 
for example, 7] is certainly as large as 4 x 10 ^ kg/W. In 
addition, we have neglected the contribution to the Birkeland 
current system driven by the pressure gradient along the 
boundary between full and empty flux tubes. This contribu- 
tion would increase the right-hand side of equation (ll) and, 
hence, lead to a larger g(0) or, in turn, a proportionately 
larger value of x. 


2. The value of A° may also be an underestimate . If, as 

Connerny et al. (1981) have estimated, we take 0 = 18°, the 

s 

value of A° is increased by a factor of thirty. 
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The distortion of Jupiter's magnetic field may also be com- 
paratively stronger than that of the earth's, which would yield a 
larger A° . One may think that 2A° 9 should, at least, be as large 
as the potential drop across the polar cap produced by field line 
merging or reconnection. Otherwise, the sunward directed convection 
resulting from such direct solar wind magnetosphere interaction 
would quench the antisunward directed planetary wind on closed field 
lines. The solar wind electric field at the orbit of Jupiter is 
400 p-V/m and if 10$ of that penetrates the magnetosphere one would 
have about 100 keV between dusk and dawn even at 9^ = 18° (assuming 
no screening out of the field by drift induced charge separation 
fields). Our estimate of A° 9 g of 1 keV may thus, indeed, be much 
too small. However, since very little is known about the occurrence 
or existence of an earth- type solar wind induced electric field, we 
will not pursue this line of reasoning. A more complete calculation 
of the equipotentials, including a possible direct solar wind 
induced dusk to dawn electric field at high latitudes will be left 
to future work. 

Finally, we mention a number of possibly observable conse- 
quences of the convection pattern discussed above (assuming a field 
of 1.5 mV/m) . 

The antisunward drift at the orbit of Io is about 700 m/ s 
increasing to 7 km/s at L 15> Such a drift may be observ- 
able either by in-situ spacecraft observations or by 


1 . 
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ground-based observations of the optical emissions from the 
Io torus . 

2. The oppositely directed flow in the cold inner torus may pro- 
duce the opposite effect, namely hotter electrons in the 
morning. However, the effect is much reduced in magnitude. 

3. Since the potential difference between dawn and dusk along 

contours of constant magnetic field is of the order of 1 MV 

at Io's orbit, one would expect a noticeable difference in 

energetic particle fluxes between the inbound and outbound 

paths of Pioneer 10, Voyager 1, and Voyager 2. Such an 

effect, corresponding to a drift contour displacement of 

about 0.1 R towards dawn, may have been observed in the 
«J 

Pioneer 10, 60 keV proton data (Randall, private communica- 
tion). Similar effects should, by the way, occur also at 
Saturn, where the physics is almost identical, even though 
the value of T| may be smaller. A search of the Pioneer 11 
Saturn data for that effect is now being undertaken. 
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FIGURE CAPTIONS 


Figure 1. A sketch of the Jovian plasma disc configuration. Flux 
tubes are emptied through a planetary wind in the tail. 

There is a boundary between full and empty flux tubes whose 
cross section in the equatorial plane is quasi -elliptical. 

The azimuthal centrifugal current flows away from this 
boundary at dusk and towards it at dawn (large arrows). 
Birkeland currents close this current system (dark arrows). 

Figure 2. Hot ion plasma densities in the front side magnetosphere 
according to Krimigis et al. (1981) and total electron densi- 
ties according to Gurnet t et al. (1981) . The boundary 
between (relatively) empty flux tubes containing low density 
hot plasma and full flux tubes containing high density colder 
plasma is indicated. 

Figure 3. Same as Figure 2 except that the data (Krimigis et al., 
1981 and Gurnett et al., 1981) are obtained in the dawn side 
magnetosphere . 
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Figure 4. A sketch of bhe equipotential contours U in the iono- 
sphere for = 4, corresponding to a flux tube content of 

, -5 , * 

6 x 10 kg/W. The contours inside of 9 = 18° are parallel 

to the dawn dusk direction. If a solar wind induced sunward 
convection (driven for example by reconvection) would exist, 
the contours would have to be appropriately rotated. 
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